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This paper investigates the possibihty of suppressing the ionization rate in a magnetostatic trap of 
metastable neon atoms by spin-polarizing the atoms. Suppression of the ionization is critical for the 
possibility of reaching Bose-Einstein condensation with such atoms. We estimate the relevant long- 
. range interactions for the system, consisting of electric quadrupole-quadrupole and dipole-induced 

dipole terms, and develop short-range potentials based on the Na2 singlet and triplet potentials. 
The auto-ionization widths of the system are also calculated. With these ingredients we calculate 
the ionization rate for spin-polarized and for spin-isotropic samples, caused by anisotropy of the 
long-range interactions. We find that spin-polarization may allow for four orders of magnitude 
suppression of the ionization rate for Ne. The results depend sensitively on a precise knowledge of 
the interaction potentials, however, pointing out the need for experimental input. The same model 
gives a suppression ratio close to unity for metastable xenon in accordance with experimental results, 
due to a much increased anisotropy in this case. 
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I. INTRODUCTION 



The experimental observation of Bose-Einstein Condensation (BEG) in 1995 by three different groups [Q-H has 
^ , fueled renewed interest in this subject, with many more groups trying to achieve the conditions under which the 
BEG-transition occurs. Stated simply, one has to reach the point where the nearest-neighbor distance between the 
^ atoms is of the same order of magnitude as their de Broglie-wavelength A. In mathematical terms, the condition to 
^ fulfill is 11 



nA^ > 2.61, (1) 



with n the atomic number density. Experimentally, this can be achieved by evaporative cooling of atoms caught in a 
trap, with typical final conditions given by a temperature of Tj, « 1 /iK and a density of n w 10^"' atoms/cm^ for 
> ' Na, and « 170 nK at n « 3 x 10^^ atoms/cm^ g for Rb. 

Both the efficiency of the cooling process and the final density achievable depend on the rates of inelastic and 
elastic collisions in the ultra-cold atom clouds used ||]. With alkali atoms, inelastic loss in binary collisions between 
' ground-state atoms is due mostly to hyperfinc-changing collisions, with rates on the order of 10~^^cm^/s for Rb [^,0. 
These rates are generallly low enough to allow condition (|l|) to be reached ||] . In addition, all the alkali atoms have 
00 ' closed optical transitions from the ground state that are accessible with commercially available lasers, allowing for 
0^ : initial trapping and cooling in a magneto-optical trap (MOT) As a result, BEG has now been achieved in Na, Rb 

and Li, while work is in progress on K and Cs. 
O , In comparison, the properties of rare gas atoms [ p^ seem less promising. As with atomic hydrogen, laser cooling 
c/j ■ and trapping starting from the ground state is not practical due to the very short wavelengths required (75 nm for 
neon). All rare gas atoms have a metastable first excited state, however, that does allow for optical manipulation. 
For He, this is the {(Is 2s) 2'^S} state, while for Ne (f^=2) through Xe (n=5) the appropriate state can be written 
{(np)^ (n-l-l)s '^P2}. Condensates of such atoms might show interesting new phenomena related to the available 
J> . electronic energy, such as collective decay or ionization. In addition, an "atom laser" |l^] based on such condensates 
may find applications that are not covered by alkali atoms: because of the large internal energy of metastable rare 
^ . , gas atoms, one might make the comparison with the high-energy photons of an XUV optical laser as opposed to those 
■ from a HeNe-laser. On the downside, this internal energy of metastable rare gas atoms is always enough to allow for 
ionization in binary collisions | [T^ . Consequently, ionization becomes a major loss process in metastable atom traps, 
with experimental rates reported to be as large as 5 x 10~^^cm'^/s for He(2'^S) |l^. Such high rates make achieving 
condition (^ impractical using the techniques currently employed for alkali atoms. 

For He, Shlyapnikov and co-workers [^,^ have shown, however, that the ionization rate can be very strongly 
suppressed by spin-polarizing the atoms. In the process 
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HeT(23S) + Het(23S) He(liS) + He+(l2S) + e 



(2) 



the initial state then has total electron spin S = 2, while in the final state spin can only be S = 0,1. Because the 
ionization process conserves electron spin, ionization is prohibited in the S = 2 (quintet) states and can proceed 
only after a magnetic interaction between the electrons causes a spin- flip during a binary encounter. Shlyapnikov 
and co-workers found that at ultra-cold temperatures the relatively weak spin-dipole or magnetic spin-spin 

interaction then limits the ionization rate in spin-polarized metastable He samples to the order of 10"^'* cm^/s |l^ , |l5| . 

For the heavier rare gases Ne through Xe, the spin-conservation rule holds also. Therefore, suppression of ionization 
by spin-polarization should be feasible as well. An added complication, however, is the fact that the unfilled (np)^ 
core now has p-character, as opposed to the He case, where it has s-character. This causes an orientation dependence 
of the electrostatic interaction potentials that is present even at large internuclear separation. During a collision, this 
causes the (np)^ cores to reorient, which, through the spin-orbit interaction, induces electron spin-flips. Since this 
interaction has an electrostatic origin rather than the magnetic one that causes spin-flips in metastable He collisions, 
one may expect it to have a much larger influence and possibly limit the ionization rate for spin-polarized samples to 
much larger values than for He. It is the purpose of this paper to estimate these effects for trapped clouds of atoms. 

Out of the various metastable rare gases, we concentrate on metastable neon because this element is being used in 
experiments presently going on in our laboratory [ p^ . In addition, the properties necessary for our calculations are 
most readily available for this species While for most alkali atoms, the properties of the atom-atom interactions 
are known quite accurately |l8[ , the situation is very different for the metastable rare gases, so that we have to make 
our own estimates. 



The {(np)^(n-(-l)s}-configuration of any of the heavy metastable rare gases contains four energy levels, with the 
•^P2-state always the lowest-lying state. Table | lists the properties of these states in neon. The four atomic fine- 
structure states give rise to a total of 144 molecular states. In what follows we will concentrate on the 25 states that 
connect asymptotically to Ne(^P2)-l- Ne(^P2), since the state is the only state that is both metastable and has 
non-zero angular momentum, thus allowing for magnetic trapping and spin-polarization. 

In the following sections we develop auto-ionization widths, short-range and long-range real potentials for the 
equivalent process of Eq. (||) for Ne('^P2). We then use them in a standard coupled-channels calculation |^ of 
the scattering matrix S. For each (positive) initial kinetic energy, we express the asymptotic solution matrix of 
the calculation as a linear combination of Bessel and Neumann functions for the radial part, multiplied by channel 
functions 4"^ denoting electronic and angular momentum states. For we have used states with well-defined total 
angular momentum quantum number P with projection M onto a laboratory axis, rotational quantum number £, 
electronic (molecular) angular momentum quantum number j and parity H, denoted as 



Since H, P and M are conserved in a collision, this has the advantage that the coupled equations can be separated 
into blocks characterized by these quantum numbers, which can then be solved separately. The 5-matrix for such a 
block will be denoted Supm- The spin-polarized and non-polarized channels can be easily identified because all the 
non-ionizing, quintet states have j = 4 exclusively. 

For any given parity, only even ^ (H = 1) or odd £{11= —1) contribute. Because the Ne("^P2)-atoms we are dealing 
with are bosons, in addition only states that are symmetric under the exchange of two atoms can contribute. This 
translates into the rule that £ + j must be even. Since parity is conserved in a collision, it follows that the non- ionizing 
quintet states {j = 4) can only couple to other even j = 0,2 states. 

Because of the presence of ionization, the 5-matrix for (in)elastic scattering is not unitary. Ionization rates follow 
from the loss of flux described by this property. For the case of a trapped cloud of atoms, the distribution of the 
direction of the initial relative velocity of two atoms is spatially isotropic. Taking this into account, we find for the 
ionization rate K^"^ in a polarized atom sample 



II. METASTABLE NEON IONIZATION RATES. 
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where the electronic angular momentum quantum number is limited to j = 4. Similarly, the rate for an unpolarized 
sample is given by 

oo 

^unpol ^ J2 1^7?°' (6) 
£=0 

oo 4 \j+£\ 

= ^EE[i + (-i)'"^]' E (2P+i){i-EEi v=(-i)w(/^'-j^)n- (7) 

i=o j=a p=\j-£\ j' i' 

A suppression ratio 77 may then be defined as 

77 = if unpoY_^pol_ 



III. LONG-RANGE POTENTIALS 

For two atoms with {(7ip)^(n+l)s}-configurations the lowest order terms in the electrostatic interaction are the 
quadrupole-quadrupolc (QQ) and dipole-induced dipolc (DD) terms |^ , ^ , with long-range behavior 

Vqq = Csi?-^ (9) 
Vdd = C&R~\ (10) 

The QQ-term is due entirely to the quadrupole moment of the tightly bound (np)^-core and is therefore relatively 
small. Its order of magnitude is given by 

a . ^yif^ (11) 

with (r|) the expectation value of the square of the orbital radius of the hole-state in the (np)^-core and q the electron 
charge. Using electron wavefunctions calculated by Hausamann ||2^ we have found (r^) ~ 1.0 a.u., leading to C5 w 1 
a.u. Numerically, we have implemented the complete state-dependence of C5 based on the theory of Refs. |po| , pT| , 
evaluated over a basis of product states | aijf'mi)| a2jf*TO2), where is the atomic electronic angular momentum 
corresponding to the atomic state labeled by a, and m its projection onto the internuclear axis. The maximum and 
minimum values of C5 that occur in the Ne('^P2)-l-Ne('^P2) subset are 0.47 a.u. and -0.37 a.u., respectively. These 
values will increase when going to the heavier rare gases Ar, Kr and Xe, in accordance with the increasing value of 



the (r^)-matrix element, which we have listed in Table III (taken from Ref. [^). According to the Table, for Xc the 
QQ-term will be about 26 times larger than for Ne. 

To obtain an estimate of the DD-term, we have used a second-order perturbation treatment of the DD-intcraction 
based on experimental values |]l^ for the {(np)^(n-f l)s} to {{n\)f{n + l)p} transition, analogous to our previous 
treatment of the DD-interaction for the {(np)^(n -I- l)s}+{(np)^(n + l)p} configuration in the rare gases ||2J|. Here 
we have included all possible product states of all {(np)^(n + l)s} and {(np)^(n -I- l)p} states, which may limit the 
accuracy of the resulting Ce values to 30-40% The orientation-dependence of the values obtained in this way 

can be understood from a semi-classical description of the DD-interaction, where the polarizability of a Ne atom is 
thought to be made up of a combination of a large, spatially isotropic value d„ « 153 a.u. due essentially to the 
(n -|- l)s valence electron, with an added small term dc ~ 1.33 a.u. determined by the polarizability of the core. Only 
this smaller term depends on the orientation of the (np)^-core, in the way described by Bussery and Aubert-Frecon 
p6[ . Such a model results in Cg values for the product state | aijY*'Tii)| Q;2j2*™2) that (for j"* — ^'2* = 2) conform to 

Ce dy^ + \dcd^{M- ml- ml), (12) 
18 

to first order in dc, which is born out by our calculated values. Since dc ~ O.Oldu for Ne, the magnitude of the 
orientation-dependent term, 4.5 a.u., is only 0.25% of the average value Cq = 1953 a.u. Again, this fraction increases 
when going to heavier rare gases, finally attaining 1.3 % for Xe (see Tables |l| and III). For He, Eq. (|l2|) is still 
appropriate for the limiting value dc = 0. 

Table |l| lists the properties of the 25 molecular eigenstates that diagonalize the combination of the DD- and QQ- 
interaction matrices with the atomic energy-matrix. From the C5 and Cg values listed, it becomes clear that, on 



3 



average, the QQ-term is quite a small effect. For one it is overwhelmed by the DD-term out to internuclear distances 
of i? « 4000 a.u. Beyond that point, its value is considerably less than 1 pK, far below the collision energy currently 
obtainable in atom traps. Therefore, the low-energy behavior of the the scattering rates will be determined by the 
DD-term. Since, however, the anisotropy due to the QQ-term (1 a.u. / R^) is of the same order of magnitude as the 
QQ-term itself, while that due to the DD-term (4.5 a.u. /R^) is only a fraction of its average value, the QQ term 
should still be expected to determine the absolute value of the ionization rates for spin-polarized atoms, which are 
governed by inelastic, spin-changing processes. 



IV. SHORT-RANGE POTENTIALS 



As mentioned in the introduction, there is little information available about the interaction potentials of the heavier 
metastable rare gases. We can obtain rough estimates of these, however, by considering the analogy with ground-state 
alkali-atoms. For the case of Ne{(2p)^3s}-fNe{(2p)^3s}, the appropriate choice is to base an estimate on the potentials 
of the system Na{(2p)^3s}-|-Na{(2p)^3s}, consisting of the Na2('^Sj) and Na2(^Sg ) states that have been studied by 
a number of authors [ pTj . 

Following Hennecart et al. ||28| , p9| , we can expand each of the four Ne{(2p)^3s} atomic states (labeled by a = 1 ... 4) 
on a basis of single-electron valence- (superscript v) and core-states (superscript c). 



aj°-^m) = X! 



1 1/2 1/2 

^ ^ ^ ^ '^m^.mj.mj (13) 

m1 = -l m5 = -l/2 mJ = -l/2 

x| (r = 1) mf)| (s- = i) ml)\ (r = 0) (m[ = 0))| (s" - i) ml). (14) 

In the following we will write | 0") as short notation for | [F — 0) (m^ — 0)) and leave out quantum numbers that do 
not change. If we write the products of the valence electron states in terms of triplet and singlet states, 

I T+i) = I +1)^ (15) 

I T_i> = I ml= ml^ -i)^ (16) 

I To) = -^1 0")J 0")2 X (I val= -1)^+ I (17) 

\S) = ^ 0")J 0")2 X (I val^ -t-i)J val= -\)- \ ml= , (18) 

then the Na2 potentials correspond to matrix elements of the atom-atom interaction Vwa given by: 

Na2(3S+) : (T, | ^Na | T^5^^, (*, J = -1,0, 1) (19) 
Na2(iE+) : (5| 1/Na I 5). (20) 

From this we can construct Ne2 potentials if we multiply VNa by a unit core-factor that is diagonal in all core quantum 
numbers. Such potentials should contain the essential features of the valence-electron exchange interactions, but do 
not describe the interaction of the core hole of one atom with the valence electron of the other. Since the radius of 
the valence electron's orbital is considerably larger than that of the core-electrons , this approximation should be 
reasonable. Also, core- valence electron terms can be added in a semi-empirical way. For an accurate description of 
the wells and turning points, however, either ab-initio potentials or experimental data would be required. 

Figure (|^) displays diagonahzed Ne(^P2)-|-Ne(^P2) potentials developed in this way, together with the Na2(^Sg ) and 
Na2(^S^) potentials on which they are based. The 9 possible, fully spin-polarized, Ne2 states that are characterized 
both by a total angular momentum j = 4 and a total electron spin S = 2, exactly coincide with the Na2('^S^) states. 
The repulsion of the valence electrons leads to a turning point i?tp ~ 8 a.u. for vanishing asymptotic kinetic energy 
in this case. 

All states conform to Hund's case (c) ||3l| , and in addition can be characterized by their total angular momentum 
j — . . . 4. Consequently, their g/u-symmetry follows immediately from j as (-l)-* . The well depth of these Ne2 
states increases with decreasing j, since more and more Na2(^S^)-character becomes mixed in. The state with the 
deepest well, of symmetry 0^, does not coincide with the Na2(^S^) state since it is not a pure singlet. In order for 
deeper wells to occur, the spin-orbit coupling of each of the Ne atoms must be allowed to compete with the Na2(^SjJ")- 
Na2('^Sj) splitting, which requires mixing with states correlating to other asymptotes than Ne('^P2)+Nc('^P2) from 
within the same configuration |(2p)^3s|-t -|(2p )^3s|. We will discuss possible consequences of including only the 



Ne(''P2)+Ne(''P2)-configuration in section VIII 
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V. AUTO-IONIZATION WIDTHS 



In scattering calculations, ionization can be introduced by adding an imaginary part to the real-valued atomic 
interaction potentials creating a complex "optical" potential 

W{R) ^V{R) -iT{R)/2, (21) 

with r(i?) the so-called auto-ionization width. 

Op de Beek et al. p3| , ^ have described a numerical method for the calculation of T{R), based on a recurrence 
recipe developed by Rico et al. jsHj. Briefly, one calculates matrix elements between initial and final states of the 
Coulomb interaction between the individual electrons involved in the ionization process. In the process considered 
here, 

Ne(3p2) + Ne(3p2) ^ Ne{^S) + Nc+{{2pf ^P)} + c" + 11.7 eV, (22) 

the dominant ionization mechanism is the so-called exchange mechanism . This means that the valence electron of 
atom 1 transfers to the available core-state of atom 2, leaving the core of atom 1 to form a positive ion. The valence 
electron of atom 2 is now no longer bound and leaves the scene. In the basis of valence-electron and core-states 
introduced in Eq. ([T^), a partial ionization amplitude 7 may be defined, given by 

7= [ dndf2mriminr-, — S ^ i ^2(^2)^^(r2)* (23) 

J 47reo \ ri - r2 \ 

with fi and r2 the position vectors of the active electrons of atom 1 and 2, respectively. The various single-electron 
valence and core states can be characterized by their orbital and spin quantum numbers with corresponding projections 
on the internuclear axis, 

vf^ (X M^' = 0; (mO^ = 0) I = i; (m,)^) (24) 
«'^cx|/2= = l;(mO^) |s^ = i;(m,)§), 

while the free electron's wavefunction "if^ is given by a Coulomb wave, characterized by orbital angular momentum 
quantum number A with projection /x, multiplied by a spin state \ sj — i;mj). Since the Coulomb interaction does 
not affect spin, Eq. ( p3| ) implies a delta-function requiring (ms)i + (jns)2 = ('^13)2 + w{. In addition, conservation of 
the azimuthal part of the orbital angular momentum requires {mi)\ + {1711)2 — {''^1)2 + 1^'^ that | |< 1. 

By using the inverse of the expansion in Eq. (^3|) on atom 1 as well as atom 2, 7 may be transformed into a 
molecular basis. From the resulting molecular ionization amplitudes 7™°'^ matrix elements of the ionization width 
r(i?) can then be constructed by summing over the available final states (represented symbolically by F), 

r{R\^^ = I ^^fITf (25) 

2Eeh'^/me F 

where i,j denote individual molecular states, rrie is the electron mass and Eg the energy carried away by the electron 
(Se = 11-7 eV for the process of Eq. (|2|)). 

We have calculated all elements of T{R) relevant for the ionization process of Eq. (|2^), using electron wavefunctions 
tabulated by Hausamann and using a modified version of the computer code developed by Op de Beek et al. 
|^,|^. Figure (||) displays the results in the same basis as used for Fig. (|l|), for internuclear distances 4 < i?(a.u.) < 12. 
As expected, the quintet states all have zero widths. Also as one would expect, the widths of the other states show 
an approximately exponentially decreasing behavior with R. At each i?, the ionization width is largest for the states 
with the largest {S — 0)-component. 

VI. CALCULATED IONIZATION RATES FOR METASTABLE NEON. 

A. Threshold behavior. 

With the ingredients developed above, ionization rates Ki for polarized and unpolarized atom clouds were calculated 
as a function of the collision energy E at ultra-cold temperatures. Figure (||) shows the result for ^ = ... 2 in the 
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range 1 /uK < E < 100 mK. The curves show the typical behavior expected from Wigner's threshold law iQ for an 
exothermic process, Kg oc E^. At very low temperature, the contribution for £ = always dominates, leading to an 
essentially constant ionization rate in this regime. As discussed in Sec. (||), due to bose-symmetry curves for odd 
partial waves are absent in the case of polarized atoms. 

The calculations of Fig. ^ show that a strong suppression, jy « 10*, is still possible notwithstanding the long- 
range anisotropy of the potentials. The resulting suppressed ionization rate, KP°^ « 5 x lO^^^cm'^/s a.t E — 1 /iK, 
would certainly allow one to reach the BEC transition, being comparable to the depolarization rate for Rb, K'^'^p"^ — 
1 X 10~^^cm^/s This result is largely due to the fact that the depolarization due to the long-range anisotropy 

is counteracted strongly at short range by the separation of the non-ionizing quintet states from the other, ionizing 
states [see Fig. (Q)]. This is illustrated in Fig. (P, which displays the rates that result from artificially replacing the 
Na2(^S^)-potential by the one for Na2('^Sj) in the calculation of the Ne('^P2)+Ne('^P2)-potentials. In this case, no 
short-range separation between ionizing and non-ionizing states occurs, with the result that ij is reduced to 77 « 170. 
This shows that the short-range behavior of the potentials is very important for the amount of suppression achievable. 

To investigate the relative importance of the quadrupole-quadrupole vs. the dipole-dipole interaction, we have 
calculated the ionization rates for polarized and non-polarized samples, with and without the QQ-term in the long- 
range potentials. Fig. (H) shows the results. From this figure, it becomes clear that the QQ-term is the dominating 
infiuence on the ionization rate of spin-polarized samples, which change by about an order of magnitude when the 
QQ-term is removed. At the same time, however, removing the quadrupole term has little infiuence on the threshold 
behavior of the rates nor on the ionization rate for unpolarized samples, confirming the dominating role of the dipole 



term on the elastic scattering already anticipated in Sec. Ill 



B. Influence of quasi-bound states. 



Since the short-range potentials for the Ne-system are not well known, we have to consider variations of the ionization 
rates due to the uncertainty in the potentials. Of particular importance are resonances due to quasi-bound states with 
total energy equal to that of the £ = quintet state. These can give rise to the formation of a relatively long-lived 
collisional complex, in which the atoms undergo a number of vibrations before they exit. In such a case, the turning 
point at small R, where ionization is most likely, may be encountered several times, leading to an enhanced ionization 
probability. To investigate the importance of such resonances, we have added a Gaussian "bump" of variable height to 
the bottom of the weU of either the Na2(^S+) or the Na2(^Sg ) potentials from which the Ne(^P2)+Ne(^P2)-potentials 
are calculated. By increasing the height of the bump, we decrease the total semi-classical phase cj) available in the 



well, effectively moving the ro- vibrational states it contains up in energy. The phase (p is calculated from |37 



/ dRj-2nV{R)/h^ (26) 



with V{R) the potential under consideration and i?tp its classical turning point for vanishing asymptotic kinetic 
energy. Figures (||) and (|^) display the variation of K^"^ and K^'^^"^ with the phase-content of the Na2(^S^) and 
Na2(^E^) potentials obtained in this way. Since we vary the phase-content by more than tt, resonances are certainly 
present in both cases. 

Figure shows that a variation of the Na2(^S^) singlet potential only has a small effect on 7^?°'. In this case, 
only the energy levels in the quintet states are unaffected by the change in phase-content, since they derive solely from 
the Na2(^i;+) triplet potential. We can conclude, therefore, that resonances in states other than the quintet states do 
not seriously alter K^"^. The most likely explanation for this is that such resonances are extremely wide, due to the 
large ionization probability in the non-quintet states, and therefore correspondingly weak. In contrast, a variation of 
the Na2(^S^) potential depth has a very large effect on K^°^, as found from Fig. (^. Here, atoms that initially have 
£ — in the quintet, are coupled to an ^ 7^ quintet state at the same total energy, which has a centrifugal barrier at 
i? ~ 75 a.u. Such a quasi-bound level may have a very sharp and therefore strong resonance, since in this state, too, 
ionization is now prohibited. During its strongly enhanced lifetime, the molecule may therefore make many vibrations 
before it tunnels through the barrier or is transferred back to the initial state. All the while, the anisotropy present in 
the electrostatic interactions weakly couples the molecule to the ionizing states, leading to a greatly enhanced K^°^. 
In our calculations, K^^]^ reaches a maximum value of 5 x 10~^"'^cm'^/s at resonance, at which point 77 « 1. Again, 
these values would be prohibitive for reaching the BEC transition. Over the whole phase-variation At/) = tt, however, 
the suppression ratio mostly takes values between 10^ and 10"*, showing that a large suppression is rather likely. 

These results clearly indicate the extreme importance of an accurate knowledge of the full potentials, or, equiv- 
alently, the scattering lengths for the process under consideration. We expect that such information will only be 
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obtained reliably from an ab-initio theoretical treatment of the short-range potentials in combination with experimen- 
tal spectroscopic determinations of the bound molecular states. Until now, this information has only been available 
for some of the alkali atoms. Since all metastable rare gas atoms other than He have more than one bosonic isotope, 
one may hope to avoid such resonances by an appropriate choice of isotope. Since there are a total of nine quintet 
states, however, with as many as 5 different potential curves if additional short-range interactions are considered, 
there may be so many different resonances that avoiding them altogether is impossible. 



C. Variation with long-range anisotropy. 



In order to make estimates of the suppression ratio for the heavier rare gases Ar thr ough Xe, we have to study 
the dependence of the ionization rates on long-range anisotropy. According to Table III , both the anisotropy in the 
dipole-dipole term as well as that in the quadrupole-quadrupole term becomes progressively larger going from Ne to 
Xe. 

A convenient way of representing the variation in C5 is by introducing a parameter AC5 = C'^°'^-C^™ equal to 
the difference of the C5 coefficients of the most and least attractive long-range potentials. Table (||) gives AC5 
= 0.84 a.u. for Ne, which may be considered its "nominal" value, which we will denote by ACs"™. We can then 
study the influence of long-range anisotropy by calculating the ionization rates as a function of ACs/ACs"™. Such 
a calculation is shown in Fig. (@), where, for clarity, we have omitted any anisotropy in the DD-term by using an 
identical, averaged Ce- coefficient for all states involved. The figure points to a quadratic dependence of ifP"' on the 
dominating anisotropy, while JC""p°i is basically unaffected by the change in C5 . This result can be explained by 
simple qualitative arguments. The anisotropic long-range interaction C5 / will cause a rotation of the molecule's 
electronic angular momentum j over an angle that is proportional to the variation AC5 in C5 , as follows from the 
equation of motion /dfi/di = r for the classical rotation angle 17 of an object with a moment of inertia / acted upon 
by a torque r (for an anisotropic potential, r is proportional to dV{R)/d^). In a binary collision of two polarized 
atoms, by analogy the initially pure quintet state will develop a projection | P) proportional to AC5 on other, ionizing 
states. Since the probability for ionization in such a collision will be proportional to (P | P), we obtain the quadratic 
dependence of ii'P"' on AC5 found in Fig. (||). 

From this figure we can therefore conclude that the limit on the suppression ratio set by long-range anisotropy 
quickly becomes stricter as the rare gas becomes heavier, because AC5 is proportional to the square of the quadrupole 
moment of the core (see Eq. (|2|)). We have listed the corresponding matrix elements in Table III for He through Xe, 
based on values taken from Ref. and indicated the relative anisotropics of all the rare gases in Fig. ||. For He, 
AC5 = since only s-electrons are involved so that no limit on rj results from the long-range anisotropy. (For He 
therefore, a first limit is set by the spin-dipole interaction For Ar through Xe, however, AC5 is appreciably 

larger than for Ne, suggesting that much smaller suppression ratios should be expected for these heavier atoms. 
Indeed, the available experimental evidence for metastable Kr |^ and Xe |3^, with, according to Table IH, approx. 
200 and 700 times larger ACs^ than Ne, respectively, showed basically no suppression at all. 



VII. COMPARISON WITH METASTABLE XENON EXPERIMENTS. 

The most extensive set of data is available for Xe, where recent experiments | p9[ | suggest that spin-polarization 
does not lead to suppression of the ionization rate for any of the available isotopes; for all bosonic isotopes, spin- 



polarization even seems to increase the ionization rate. From the discussion in Subsection VI C one would already 
expect 77 to be less than 10 in this case. However, additional dynamical effects due to the larger mass of Xe, and 
the difference in short-range potentials compared to Ne should also be taken into account. These changes can be 
easily incorporated into our model by substituting Cs potentials |4^] for the Na potentials used for Ne, and using the 
appropriate mass. We have used the same auto-ionization widths as for Ne, however, since auto-ionization widths for 
Xe are not available. 

Figure ^ shows the result of such a treatment for a collisions energy of 1/iK. Again we varied the amount of phase 
available in the triplet-potential to obtain the variation of the ionization rates for polarized and unpolarized clouds 
of (Xe) atoms. Indeed, suppression by a factor of -q «10 is still possible, but over 60% of the 0-range in Fig. ||, 
?7 < 1. In fact, the 0-averaged value is now {rf) = 0.7, showing that for Xe ionization suppression by spin-polarization 
is very unlikely. Here, the QQ-term mixes up the polarized and unpolarized states so thorougly at long-range that 
they basically ionize at the same rate. The fact that 77 < 1 can then be considered a quantum-statistical effect: for 
polarized atoms, the ^ = partial wave necessarily represents a greater fiux than for unpolarized atoms, since in the 
former case the ^ = \ partial wave is absent. 
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VIII. DISCUSSION. 



In this paper we have developed a semi-quantitative description of ionization in metastable neon collisions at very 
low temperatures. In fact, this description holds for all the heavy metastable rare gas atoms, but we have concentrated 
on neon because of ongoing experiments in our laboratory as well as the availability of the properties needed as input 
to the model. We have put special emphasis on obtaining a lower limit on the ionization rate for spin-polarized 
atoms, in order to investigate the possibilites of reaching Bose-Einstein Condensation with metastable rare gas atoms. 
While we are not able to provide a definitive answer, we do provide a framework of potentials and auto-ionization 
widths that allow one to investigate the system in a qualitative way. It has become clear that such a limit depends 
strongly on the occurrence of quasi-bound resonances in the potentials for fully polarized atoms with £ > 0. While 
it is possible to obtain reasonable estimates of the anisotropy of the long-range part of these potentials, sufhciently 
accurate information on the short-range parts is not available at this moment. Such potentials may become available 
in the near future, however |41|. Our calculations show that, under favorable conditions, the long-range anisotropy 
still allows for very effective suppression (up to 10* times) of ionization by spin-polarization, at least for metastable 
Ne. This suppression becomes less effective as the atom becomes heavier, with a limit of 77 « 10 for Xe. 

Additional couplings may exist between the quintet states and the ionizing states that we have not considered. 
Such couplings may occur at smaller R, where core-valence interactions become more important. The interactions 
that will have an effect on the ionization rate for spin-polarized atoms must create a direct coupling between a quintet 
state and other, ionizing states. For Xe, e.g., one should expect a large effect of the second-order spin-orbit splitting 
from the experience with Cs [^,0. A simple splitting of the various quintet potentials, however, due to a dependence 
of the potentials on a cr or tt orientation of the cores, does not have such an effect, at least not for small magnetic 
confining fields. This is because metastable rare gas atoms in general lack hyperfine-structure, so that exothermic 
exit channels do not exist. As a result, depolarization through such a mechanism is prohibited. We have checked this 
by adding a phenomenological a/n core-splitting to our potentials, with an order of magnitude comparable to that 
found in the Ne+Ne^-system E3]. The ionization rate for polarized atoms was found to increase only slightly with 
the magnitude of this splitting. 

A large magnetic confining field, however, could cause Zeeman splittings greater than the centrifugal barrier present 
in the £ > quintet states. This would open up new channels for depolarizing collisions that require only a splitting 
within the quintet states. Given the height of the centrifugal barrier, 5.8 mK for Ne and £ — 2, such processes should 
be expected to be fully developed at magnetic fields of approximately 60 Gauss. It may be profitable, therefore, to 
use the minimum practical trapping field for metastable rare gas atom traps. 

Another factor that we have not taken into account is that curve-crossings of the Ne('^P2)-l-Ne('^P2) states can occur 
with other molecular states from the Ne{(np)^(n-|-l)s}-|-Ne{(np)^(n-|-l)s} system at some intermediate range. For 



Ne, it is possible to develop potentials for this full system with the recipe developed in Sec. IV. The result of such 
a treatment is shown in Fig. jl^). In this case, the Ne('^P2)-l-Ne('^P2) quintet states indeed cross several states with 
other asymptotic limits. Closer inspection has taught us that some of these are indeed avoided crossings, with the 
coupling due to the DD- and QQ-interactions. At the crossing points, a second order treatment of the DD-interaction 
in fact no longer suffices. Flux entering on one of the quintet states can now be diverted onto curves that allow 
ionization to occur at short range. In order for the effect of such crossings to be taken into account, an i?-dependent 
treatment of the DD- and QQ-interactions would be necessary, and in addition a much larger system of coupled 
equations would have to be solved. We consequently leave such a treatment to a future paper. 
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TABLE I. Spectroscopic notation, equivalent Russel- Saunders LSJ notation where applicable, atomic energy and electronic 
angular momentum quantum number J for the four states of Ne in the {(2p)^3s}-configuration 



Spectroscopic 


LSJ 


Energy 


J 


notation 


notation 


(cm-^) 




3s' [1/2] 1 




135888.7137 


1 


3s' [l/2]o 




134818.6405 





3s [3/2] 1 




134459.2871 


1 


3s [3/2] 2 


^P2 


134041.8400 


2 



TABLE IL Hund's case (c) symmetries, Cgand C5 coefficients of the 25 adiabatic molecular states asymptotically connected 
to Ne(^P2)+Ne(^P2), derived from the Na2(^E+)and Na2(^E+)-potentials with the formalism of Sec. The de generacy of 
all f2 = states is 1, that of all others is 2. 



n 


C5 (an) 


Ce (au) 


3 


-0.3700 


-1953.24 


2 


-0.2878 


-1951.33 


2 


-0.2673 


-1952.60 


1 


-0.2028 


-1951.23 


0" 


-0.1525 


-1951.64 


1 


-0.1318 


-1951.13 


0+ 


-0.0679 


-1954.13 


3 


0.1233 


-1953.36 


2 


0.1233 


-1952.73 


1 


0.1617 


-1952.78 


1 


0.1730 


-1953.52 


0+ 


0.2173 


-1955.58 


4 


0.2467 


-1955.54 


0" 


0.3992 


-1951.48 


0+ 


0.4673 


-1952.38 



TABLE in. Atomic and core polarizability, and dc, respectively, and expectation value of the square of the orbital radius 
of the hole-state in the (np)^-core, (taken from Ref. [Q) for metastable rare gas atoms Ne through Xe. The atomic polarizability 
corresponds to the value listed by Ref. [^ for configuration {(np)^(n-f l)s}, while the core polarizability corresponds to the 
ionic configuration (np)^. Since dc/dv <C 1 always, one may equate the polarizability of the (n -I- l)s-valence electron alone 
with the atomic polarizability. The last two columns give the anisotropy parameters ACe and AC5 relative to those for 
Ne('^P2)(AC5 — 0.84 a.u., ACe ~ 4.5 a.u.). For clarity, the table also lists the anisotropy parameters for He as being zero: in 
this case the core has a spatially isotropic (Is) configuration. 



rare gas 


principal quantum 


dv 


dc 




anisotropy ACe 


anisotropy AC5 


atom 


number n 


(a.u.) 


(a.u.) 


(a.u.) 


relative to Ne(^P2) 


relative to Ne(^P2) 


He 
















Ne 


2 


153 


1.33 


1.23 


1.0 


1.0 


Ar 


3 


305 


7.90 


3.31 


3.0 


7.3 


Kr 


4 


363 


12.8 


AA6 


4.1 


13.2 


Xe 


5 


511 


23.1 


6.28 


5.2 


26.1 
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FIG. 2. Ab-initio auto-ionization widths F for the Ne( P2)+Ne( P2)-system as a function of internuclear distance R. The 



widths are given in the same, adiabatic basis as used for Fig. 
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The inset shows the widths on a logarithmic scale. Note the 
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FIG. 3. Partial ionization rates for trapped Ne(^P2)-atoms as a function of the collision energy E. Rates are given for 
unpolarized {K^"^°'') as well as for polarized atoms {K^°^) for angular momentum quantum numbers I = Q . . .2. The rates 
show an obvious Wigner-law behavior. The rates for polarized atoms are about 10* x lower than those for unpolarized atoms. 



C/3 



-I— > 

o 

• 1— I 
-I— > 

N 

• 1— I 

o 



10- 


10 


10- 


11 


10- 


12 


10' 


13 


10' 


14 


10' 


15 


10' 


16 


10' 


17 



: 1=0 ... — r--:::... ^ 

•• •••• *•• 

.* / \ •••• 


1 1 1 M IM| 1 1 M 1 Ml| = 

un- 1 
polarized 


?/=i./ / / \ 




r ,: ^..^^^^i^-'^N^^^^^^^^v^ 


polarized 


[ / /i=o 




U=2/ / 

1 1 III 1 1- 1 1 1 1 1 1 


1 1 



10"^ 10-2 10-^ 10^ 10^ 10^ 10^ 10"^ 



kinetic energy E (mK) 

FIG. 4. Partial ionization rates for trapped Ne(^P2)-atoms as a function of the collision energy E. In evaluating the 
potentials for the system, the Na2(^E^)-potential was artificially replaced by the Na2(^E^) potential, removing the large 
energy-splitting between ionizing and non-ionizing states. The resulting suppression ratio is much less than for the case of Fig. 
(i). 
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FIG. 5. Partial ionization rates for trapped Ne(^P2)-atoms as a function of the collision energy E, calculated including 
or excluding the effect of the quadrupole-quadrupole (QQ) interaction on the potentials. For the polarized case the results 
are labeled "DD + QQ" and "DD only", respectively. In general, these rates drop by about an order of magnitude when the 
QQ-term is omitted. For the unpolarized case, the solid lines give the result including, and the triangles the result excluding 
the QQ-term. Little difference is observed in this case. 
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FIG. 6. Partial ionization rates for trapped Ne(^P2)-atoms as a function of the depth of the Na2(^E^)-potential from which 
the potentials of the Ne(^P2)-|-Ne(^P2)-system are calculated. The depth of the Na2(i^Sj^)-potential is given in terms of the 
semi-classical phase [see Eq. (p^)]. Only relatively small variations in the rates are observed, indicating the minor influence of 
quasi-bound resonances in the ionizing states on the ionization rate of the non-ionizing, polarized quintet states. The collision 
energy is i5 = IfiK. 



13 



^ 10' 



s 

u 

i4 

B 

a 

N 

'S 

o 



10 



12 



10" 



10 



10 



10 



18 



■21 



■24 



non-polarized 



polarized 




14.0 14.5 15.0 

triplet phase (7t radians) 



14.0 14.5 15.0 

triplet phase (n radians) 



FIG. 7. Partial ionization rates for trapped Ne(^P2)-atoms as a function of the depth of the Na2(^Sj )-potential from which 
the potentials of the Ne('^P2)+Ne(^P2)-system are calculated. The depth of the Na2(^S+)-potential is given in terms of the 
semi-classical phase (j> (see Eq. (p6|)). Very large variations in the rates are observed, indicating narrow quasi-bound resonances 
occurring in the non-ionizing quintet states. These may adversely affect the possiblities of achieving BEG with metastable rare 
gases. The collision energy is i5 = IfiK. 




FIG. 8. Partial ionization rates for trapped Ne(^P2)-atoms as a function of the normalized long-range anisotropy parameter 
ACs/ACg^"™. A quadratic dependence of the rate for polarized atoms on the anisotropy parameter is observed. Arrows 
indicating the estimated anisotropy parameter for metastable He, Ne through Xe (listed in Table III), illustrate the decrease 
in ionization suppression ratio with increasing rare gas atom mass. The collision energy is _E = 1/jK. 
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FIG. 9. Same as Fig. M but for Xe instead of Ne, replacing the Ne mass by that of Xe, and the Na-potentials by those of 
Cs. The larger anisotropy in the case of Xe leads to much greater polarized ionization rates than for Ne, often exceeding those 
for unpolarized atoms. A suppression ry = 10 is possible, but values rj < 1 are more likely. The collision energy is iS = l/iK. 
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FIG. 10. Adiabatic potentials for the Ne{(2p)^3s}+Ne{(2p)^3s}-system including all fine-structure states. The zero-point 
of the energy scale corresponds to the asymptotic energy of the Ne(^P2)+Ne('^P2)-configuration. Compared to the re- 
duced set of Fig. (|l|), deeper wells (resembling more closely the Na2(^E^) state) have developed in the states connected 
asymptotically to Ne(''P2)+Ne(''P2), indicating the relaxation of the asymptotic spin-orbit configuration. Crossings of the 
Ne(^P2)-|-Ne('^P2)-system with other asymptotic configurations are indicated as well. Such crossings may lead to reduced 
ionization suppression by funneling flux onto ionizing potentials. 
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